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We study the cosmological evolution of extended branes in 6D warped flux compactification mod- 
els. The branes are endowed with the three ordinary spatial dimensions, which are assumed to be 
homogeneous and isotropic, as well as an internal extra dimension compactified on a circle. We 
embed these codimension f branes in a static bulk 6D spacetime, whose geometry is a solution 
of 6D Einstein-Maxwell or Einstein-Maxwell-dilaton theories, corresponding to a warped flux com- 
pactification. The brane matter consists of a complex scalar field which is coupled to the bulk U(l) 
gauge field. In both models, we show that there is critical point which the brane cannot cross as it 
moves in the bulk. We study the cosmological behaviour, especially when the brane approaches this 
critical point or one of the two conical singularities. In the present setup where the bulk geometry 
is fixed, we find that the brane cosmology does not coincide with the standard one in the low energy 
limit. 

PACS numbers: 04.50.+h, 98.80.Cq 



I. INTRODUCTION 

Braneworld models have been studied actively in the 
last decade (see e.g. jl| for a few reviews on this topic). 
Among models which take explicitly into account the 
self-gravity of the branes, i.e. the backreaction on the 
geometry due to the presence of the brane, most efforts 
have been devoted to the study of gravity and cosmology 
in five dimensional (codimension 1) braneworld models, 
before (see @) and especially after the stimulating pro- 
posals by Randall and Sundrum [3J- However, if one is 
interested in a braneworld model in a six or higher di- 
mensional bulk spacetime, one needs to consider a brane 
with codimension higher than 1. 

Among higher codimensional braneworld models, some 
significant attention has been paid to codimension 2 
branes (see e.g. Ref. [H and references therein), because 
they possess the property that, even with some non zero 
vacuum energy, the brane geometry can be flat, while the 
bulk geometry surrounding the brane is characterized by 
a deficit angle. Codimension 2 branes have thus been 
studied in the context of the cosmological constant prob- 
lem [1, H, 0| , but they are also interesting on their own 
as the simplest case beyond codimension 1 branes. 

One difficulty however with codimension 2 brane is 
that the energy momentum tensor is in general restricted 
to be that of a pure tension. To consider more general 
matter, such as would be relevant to study cosmology in 
the brane, one needs to regularize the brane. One well- 
motivated way of regularization is to replace the codi- 
mension 2 brane by a thin ring-like codimension 1 brane 
wrapped around the symmetry axis, as was firstly sug- 



gested in Q in the context of the "rugby-ball" braneworld 
model. The rugby-ball model has a 2D compact bulk 
compactified by a magnetic flux and codimension two 
branes are located at poles. This model was originally 
introduced in order to resolve the cosmological constant 
problem [f| , but it was also pointed out that the so-called 
selftuning property did not actually work [f| 0] (see also 
Q for a review). 

In the present work, we consider static warped 6D 
braneworld solutions [l(J, EH with a U(l) gauge field cou- 
pled to gravity, which are generalizations of the rugby- 
ball model. Only static regularized branes have been 
considered in these geometries so far as now [lU, E3, EH ■ 
The goal of the present work is to study the cosmologi- 
cal behavior of branes in such 6D braneworld models, i.e. 
to consider the motion of the branes in these given bulk 
geometries. 1 We basically follow the well-established 
method to discuss cosmology in 5D braneworld models, 
using the static bulk point of view [HI, El, ES EH rather 
than the brane-based approach 0] . In order to satisfy all 
the junction conditions, we assume, like in Q that the 
branes contain a complex scalar field but we allow for a 
time evolution of its radial component. 

The paper is organized as follows. In Sec. II, we give 
a general formalism to discuss cosmology on the brane. 
In Sec. Ill, we focus on the cosmological evolution of the 
brane in the models based on the 6D Einstein-Maxwell 
theory pjj. In Sec. IV, we discuss the case of the 
6D Einstein-Maxwell-dilaton model (more precisely the 
bosonic part of Nishino-Sezgin supergravity [HI]) [ll[ as 
the same manner as the line of the previous section. In 
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1 In Ref. 0, cosmology in such models has been discussed by 
employing a different way of regularization. 
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Sec. V, we give summary and conclusions before closing 
the article. 



II. BASIC FORMULATION 

Before considering explicit models, we study in this 
section the general problem of the motion of a codimen- 
sion 1 brane embedded in a 6D static geometry with axial 
symmetry in its two internal dimensions. We start from 
a static 6D metric of the form 

ds 2 = A 2 {w)(-dt 2 + dx 2 ) + B 2 {w)dw 2 + C 2 (w)d9 2 (2A) 

Here 9 is implicitly assumed to be an angular coordinate, 
associated with a compactificd internal dimension. A(w) 
is assumed to be a monotonically increasing function of 
w. We embed into this static bulk a codimension 1 brane, 
which we assume to be homogeneous and isotropic along 
the three ordinary spatial dimensions as well as along the 
angular coordinate 9. The position of the brane is thus 
simply characterized by its coordinates t(r) and w(t), 
given for example as functions of the proper time r. By 
definition of the proper time, one has the following rela- 
tion 



A 2 i 2 + B 2 w 2 



■1 . 



(2.2) 



where a dot denotes a derivative with respect to r. The 
above relation comes from the normalization of the brane 
5-velocity, whose components are given by 



on the left hand side denote the jump across the brane 
of the extrinsic curvature tensor, that is the difference 
between its values on the two sides of the brane. 

In our case, since we assume homogeneity of the brane 
matter as well as isotropy in the three ordinary spatial 
dimensions, the brane energy momentum tensor is nec- 
essarily of the form: 



T ab = £ u a u b + p hah 



(2.8) 



where h ab is the three-dimensional metric in the ordinary 
spatial directions and 9 a is the unit vector along the in- 
ternal 9 direction. E is the energy density and p and pe 
are respectively the ordinary and internal pressures. 

One consequence of the brane junction conditions is 
the generalized energy conservation law on the brane, 
which can be expressed as 

E + 3- (E + p) + ^ (E + Pe ) = [T ab n a u b ] , 

,+ - (^ b n a u b ) (2.9) 



[T ab n a u b ] = (T ab n a u b ) 



where Tab denotes the bulk energy-momentum tensor, 
and all quantities are evaluated at the brane location. 
The above expression differs from the usual cosmologi- 
cal conservation law in two respects: first, there is an 
additional term due the evolution of the internal spatial 
direction; second, one finds in general on the right hand 
side a brane-bulk energy exchange term, similar to what 
can been found in 5D brane cosmology (see e.g. |13|). 



u t = t u w =w u t = ~A 2 t u w = B*w . (2.3) 
The unit normal vector to the brane is defined by 



n w = -i n* = ^rw 
B A 

n w = ABt n t = —ABw. 



(2.4) 



If q a b denotes the induced metric on the brane, the ex- 
trinsic curvature tensor is defined by 



K ab = q a c V c n b = ^C n q ab , 



(2.5) 



where £ n is the Lie derivative with respect to the normal 
vector n a . In order to write explicitly the junction condi- 
tions, it is more convenient to introduce the combination 



Kq ab , 



(2.6) 



where K = q ab K ab . The junction conditions are then 
given by 



K„ 



1 

Mi 



(2.7) 



where T ab is the energy momentum-tensor confined on 
the brane and Mq is the 6D Planck mass. We shall set 
Mq = 1 unless stated otherwise explicitly. The brackets 



III. EINSTEIN-MAXWELL MODEL 
A. Bulk spacetime in 6D Einstein-Maxwell theory 

We now discuss the dynamics of an extended brane in 
the 6D Einstein-Maxwell theory. The action in the bulk 
is given by 



d b x^ 



i 



:F AB F 



AB 



Ao) , (3.1) 



where gAB is the 6D metric, Fab — 28[aAb] is the field 
strength associated to the U(l) gauge vector As, and A 
is a cosmological constant. 

By a double Wick rotation of the 6D Reissner- 
Nordstroem solution, one can obtain the solution p^ |: 



ds 2 



9 rjnndx^dx 1 



dp 2 
Fjpj 



c 2 F{p)d V 2 



F, 



where 



bo 



\2c 2 p' 



2n6 



n 3 



(3.2) 



(3.3) 
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Assuming that F has two positive roots p_ and p + such 
that < p_ < p + , one can replace the two parameters 
60 and po by p+ and a = p-/ p+. Using F(p±) — 0, feo 
and po are given by 



geometry at the two poles w = 1 and w = a is now 
smooth, i.e. with no deficit angle. This implies [HI, [l3j] 
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Mo 



10 P+ 1 



or 



(3.4) 



20(l-a)(l-a 3 ) 
5 - 8a 3 + 3a 8 ' 
20(l-a)(l-a 3 )a 4 



(3.10) 



It is also convenient to replace the coordinates p and ip 
by the new coordinates w and 8, defined respectively by 



p = p+w. 



<p = 



A oP+(l - a) 
so that the bulk metric now reads 
Mw 2 



(3.5) 



5a 8 - 8a 5 + 3 

Moreover, the continuity of the (0, 9) component requires 

t := Rqcq = R+c+ = R-C- . (3.11) 

For a schematic picture of the bulk configuration, see e.g., 
Fig. 1 of Ref. 

The gauge field in the bulk is solved as 

2S£ 



dsi 



R 



(l-a) 2 /H 



4f{w)d9< 



+ p\w 2 -r] ilv dx ll dx v , 



(3.6) 



2S£ 1 1 



3(1 — a) \w 3 



1 



Q" 



(3.12) 



with 



1 



10 



1 a 3 (l - a 5 ) 



1 1 



5(1 -a) 2 

and the field strength is now given by 



1-a 3 



ui 3 1 



(3.7) 



expression valid respectively in the northern and south- 
ern regions. Note that in the region of the "0"-bulk the 
bulk gauge fields are doubly defined. This is consistent 
only if the difference is simply a pure gauge term, which 
implies the Dirac quantization condition 



N 



2eS£(l - a 3 ) 
3(1 -a)a 3 



N = 0,±1,±2,- 



(3.13) 



with R\ = 



2c RqS 
(1 — a)w 4 

(2A )- X . 



S = S{a) 



/3a 3 (l 



a ,J 



5(1 -a 3 ) 



(3.8) 



The singularities p = p_ and p = p + now correspond 
to w — a and w = 1 in the new coordinate system. By 
expanding the above metric around these two singular- 
ities, one finds that they are in general conical singu- 
larities with some corresponding deficit angles. These 
conical singularities can be related to the presence of two 
codimension-2 branes. As first suggested in [H, one way 
to regularize these codimension 2 branes is to replace 
them by two codimension 1 branes, with one spatial di- 
mension compactified on a circle at the respective posi- 
tions W- > a and w + < 1. These codimension 1 branes 
are the boundaries of regular spherical caps, which end 
smoothly the spacetime. The spherical caps can be de- 
scribed by solutions similar to that of the main bulk, i.e. 



B. Junction conditions 

The matter content of the brane must be compatible 
with the bulk solution, i.e. must satisfy the junction 
conditions for the metric and the gauge field. In our 
case, we consider a complex scalar field coupled to the 
bulk gauge field with the following action 



Sb,± = - 



rf 5 .; 



^q{V(\4>±\) 

/ =gfv(«±)4 



eA a )(d a a ± -eA a ) 



(3.14) 



where 



d s 1,(i) 



t w9 - 



o2 [ 4d W 2 

U(l-a)3/(«0 
p\w i f]ii,vdx il dx v , 

2gRiS 
~ (1 - a)w i 



cjf(w)d9 2 



(3.9) 



± = $±(r)e I<T± ^- 1 . This is a natural extension of 
the brane action adopted in [1, 0j|, [U[ for static branes, 
with the difference that which was frozen at a fixed 
value in the static case, must now be promoted to a time- 
dependent field. Hereafter, we will often omit the brane 
subscripts ± for brevity. The above action implies that 
the energy density, the external and internal pressures, 
defined earlier in (|2.8[) . are given explicitly by the follow- 
ing expressions: 



where the index / takes the values I = + and / = — to 
denote respectively the northern cap, w+ < w < 1, and 
the southern cap, a < w < w_. Note that each bulk 
region is endowed with different cosmological constants 
A±. The parameters c + and c_ are chosen so that the 



Po 



-y($) 



-<S> 2 C- 2 {d e a-eA e ) 2 + -$ 2 , (3.15) 



<S> 2 C- 2 (d g a-eA g ) 2 
<S> 2 C- 2 (d e a-eA g ) 2 



$ 2 (3.16) 
$ 2 .(3.17) 
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To discuss the dynamics of the (±)-branes, we can now 
apply the general formalism presented in the previous 
section. Comparing the metric ()3.6() with Eq. (|2.1j) . it is 
easy to identify 



A(w) 2 =p 2 + w 2 , B(w) 2 
C(w) 2 = £ 2 f(w). 



4R 2 



(l-a) 2 /H 



(3.18) 



Note that the warp factor A is a monotonically increas- 
ing function with respect to w as assumed in the previ- 
ous section. Since we are interested in the cosmological 
evolution of the branes, it is convenient to reexpress the 
various bulk components as functions of the brane scale 
factor 



o(t) = A(w b (r)), 



(3.19) 



instead of the radial coordinate w. In the following, we 
will use the function C (a) defined by substituting in C (w) 
the expression of w as a function of a, as well as 



\dw) (i-«Kv7' 



(3.20) 



where the right hand side is evaluated at the brane posi- 
tion. 

Defining the jump across the branes as [T]± :— ±(T± — 
T ), the junction conditions for the metric [K a (,] = —T ab 
are given by 



/3 C 







■ 1 " 


X 


+ (« + fa 2 ) 








X. 



X 



B' . 
B 



r 1 

x 



P9 



/3 C'\ 
L + c) 



X 



= p, (3.21) 
(3.22) 
(3.23) 



where X := y B~ 2 + a 2 and the prime ' denotes the 
derivative with respect to a. These relations generalize 
the expressions given in [![ for a static brane. Note that 
in the static case, the left hand sides of (|3.21|) and (|3.23|) 
are identical, which implies p — — S. Note also that in the 
cosmological case, the topological condition is the same 
as in the static case (see Appendix A-l for details). 

The last junction condition (|3.23[) implies the general- 
ized Friedmann equation 

a 2 p 2 



16tt 2 £ 2 /(3+ Ija) 2 



n 2 £ 2 f 



2j a ) (B± 2 ~ Bq 2 ) 2 



B\ 



B n 



P 2 



2a 2 



(3.24) 



where p := 2Tr£y /r fT,, obtained by integrating S over the 
internal dimension, is the 4D effective energy density. 
On the right hand side, one finds a term quadratic in 



the energy density, which is familiar in brane cosmology 

0, as well as a term which goes like l/p 2 and which is 
characteristic of non ^-symmetric brane cosmology [l8| , 

1. e. when the two bulk regions surrounding the brane 
have a different geometry. 

It is instructive to study the low energy behaviour of 
the Friedmann equation by considering an expansion of 
Eq. (|3.24p around a static brane configuration at a = 
ao such that d| a=Qo = 0. From Eq. (|3.23[) . the energy 
density at a — clq is given by 



S 



for the (+) brane and 
for the (— )-brane. For the (-h)-brane, we obtain 



a— <1q 



So 



9L 
~c 



(3.25) 



(3.26) 



off r ~ 



Sp + 0(Sp 



7,2 \ 



a 2 3 

where the effective gravitational coupling is given by 
3(1 - a)p + 



(3.27) 



87rG c ff — 



2ir(R 



R+)£M£a(3 + ^f) 



(3.28) 



We observe that the gravitational coupling depends 
strongly on the scale factor and is even negative for 

( 3(1 — a 8 ) \ 

ao > p+ ( 8(i-a3) J ' name ly m t ne vicinity of the 
northern pole. Therefore, one cannot expect the cosmol- 
ogy governed by the Friedmann equation (|3.24|) to co- 
incide with standard cosmology in the low energy limit, 
in contrast to the 5D Randall-Sundrum cosmology. The 
conclusion here does not depend on the choice of brane 
matter. A similar expression for the effective gravita- 
tional coupling for the (— )- brane is obtained by replac- 
ing i? and R + in Eq. (|3.28[) with, respectively, R- and 

Rq. 

Now, we consider more specifically the cosmological 
evolution when the brane matter consists of a complex 
scalar field. We must also take into account the junction 
condition for the gauge field which is given by 



J"(o) X 



-e<£ 2 (d e <T - eA g ) , 
where we have defined the function T(a) 

2c I R I S{a) fp , x i 



(3.29) 



T{a) 



dA 
dw 



F„ 



I A=a 



(1 - a)p 



V fP±Y 
_l_ V a / 



(3.30) 

Let us now combine the four junction conditions above. 
By subtracting Eq (|3.2ip and Eq. (|3.22| ), and then using 
Eqs. (|3.16p . (|3.17[) and (|3.29p . one gets the expression 



X 



-e 2 <f> 2 



f 2 \C 



\C a) 



(3.31) 
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where the right-hand side depends only on a and <f> but 
not their derivatives. By substituting this expression 
back in the junction condition (|3.23p . one finds that the 
energy density of the brane £ is given explicitly in terms 
of a and <& only: 



-V— -- 

T 2 \C' + a)\C ~ a 
e 2 $ 2 (l- a) 2 a s f(a) //'(a) 



4S 12 



P 6 + 



2/(a) 



;)( 



/» 3 
2/(o) a. 

(3.32) 



C 2 (7 



(3.33) 



This implies that the part of the energy density that de- 
pends only on the scalar field $ is constrained to be a 
specific function of a and $, namely 

(3.34) 

By subtracting Eq. (|3.2ip from Eq. (|3.23p . we also obtain 



o H — —a' 
B 



(3.35) 



which governs the brane acceleration. 

The bulk a < w < 1 can be divided into two regions, 
where the branes will be confined: a brane in a given 
region cannot move to the other one. This can be seen 
by noting that the sign of 



[X}±=±{ 



V+/(a) 



'(1 



4i? 2 



(3.36) 

is determined by the relative values of Ri on the two sides 
of the brane and thus cannot change. From Eq. (|3.3ip . 
one sees that the sign of [X] is negative for a < w < w\ 
and positive for w\ < w < 1, where 



Wi 



8(1 -a 5 ) 
5(1 -a 8 ) 



1/3 



(3.37) 



QL 

c 



which is plotted in Fig. 1, 
W\ represents a limit that 



-5 - 
-10 - 

-15 - 
-20 - 
-25 ■ 



One can even get a more precise constraint on the energy 
density in the brane by noting that the substitution of 
(|3.3ip into the gauge field junction condition (|3.29|) yields 



FIG. 1: Numerical plot of C'/C — 1/a in the Einstein- 
Maxwell model is shown as a function of a/p+ for a — 0.1. 



changes. Since the sign of p — pg cannot change, we see 
that the brane cannot cross w\ because the junction con- 
ditions cannot be satisfied on the other side. By noting 
that i?+ < R-, which follows from the relations (|3.11|) 
and (|3.10[) . one finds only two possibilities for the loca- 
tions of the branes. The first possibility is to put the 
(+)-brane in the northern region, i.e. w + > wi, and the 
(— )-brane in the southern region, i.e. W- < w\. Another 
possibility is to put both branes in the nothern region, 
i.e. w+ > W- > u>\. It is however impossible to put both 
branes in the southern region. 

Finally, it is worth noticing that the brane energy den- 
sity S is negative when w\ < w < u>2, where 



w 2 := 



3(1 -a 8 ) 
8(1 - a 3 ) 



1/5 



is the value at which the function ^7- + - vanishes. Note 
that, whereas the coefficient in front of p 2 in the Fried- 
mann equation Eq. (|3.24p as well as the effective grav- 
itational coupling diverge when one approaches W2, the 

energy density Eq. (|3.32p tends to zero as well so that 

±he overall term does not diverge: the brane can cross 
a- 2 smoothly W2 , simply with a change of sign for the energy 
density. 



(3.38) 



is the value at which 
vanishes. Consequently, w 
the branes cannot cross. 

The existence of this critical value is a consequence of 
the junction conditions and our choice of brane matter. 
Indeed, on the one hand, we always have p > pg with 
the complex scalar field that was assumed. On the other 
hand, the extrinsic curvature tensor is dictated by the ge- 
ometry and wi is the value where the sign of K\ /3 — Kg 



C. Brane dynamics 

To summarize, we have shown in the previous subsec- 
tion that, after eliminating the explicit dependence on 
the brane Goldstone mode a and vector potential Ag via 
the Maxwell junction condition, the other three junction 
conditions yield 



£0 



a 2 - JB^ 2 + a 2 = $ 2 G(a) , 



$ 2 V2 



jVo $ 2 + n$)) =K(a) 



(3.39) 
(3.40) 



/- d a f. 2 \f 1 



B2 2 + a 2 



B\ 



= $ 2 . 



(3.41) 
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where 



K(a) 



G(a) 



8S 2 p 
e 2 (l- 



2/ 



AS 2 (A 



fwa I — - - - 

V 2/ a 



(3.42) 



The above equations apply to both branes: for the north- 
ern (+)-brane, the subscripts "1" and "2" correspond to 
the "0" and " +" bulk solutions while for the (— )-brane, 
they correspond to the " — " and "0" bulk solutions, re- 
spectively. 

The first equation p.39[) implies 



D. Behaviour near the "barrier" 

As we saw earlier, w — w\ represents an uncrossable 
barrier for the branes. It is thus interesting to study the 
cosmological behaviour of a brane that approaches this 
critical value. In the limit a — > a\ := w\p + , the function 
G vanishes. Keeping in the differential equation (|3.46p 
the terms that dominate in this limit, one obtains 



\~da~) 



8$ 6 G 3 



{B 2 



■ 1 d$ 
. $ da 



d a G d a f 



G 



2/)}' 
(3.48) 



1 



4$4 G 2 
_ -2 _ 



$ 4 G 2 - (Sf 
W(a,$) 2 . 



Bo 



3> 4 G 2 - (Sf 1 - B 2 



where one can substitute G(a) = G'(a\)(a — a\). 
One can then easily find the approximate solution 



(3.43) 



$ ~ $i-|-$i(a-oi) 



G ,d 



By noting that 



{b^-b^Y 



n 2 



$ 4 G 2 + b: 



b; 2 



Bo 



n 2 = 



$ 4 G 2 



2$ 2 G 
Bo 



Br 



2$ 2 G 



(3.44) 



\a— a± 

(3.49) 

The potential Eq. (|3.47p approaches the constant value 

(3.50) 



V 



i$ 2 G'( ai ). 



and that the acceleration of the scale factor can be reex- 
pressed as 



d = -a = a—H = n(— + — — 
dt da \ da <9<I> da 



(3.45) 



we find, after some manipulations, that (I3.41| . divided 
by a 2 , can be reduced to 



(-) 

\daJ 



4$6 G 3 



where we have used Eq. (|3.49|) . The above analysis ap- 
plies for the (+)-brane (and possibly for the (— )-brane) in 
the northern region when the brane approaches a\ from 
above. It is also valid in the southern region when the 
(— ) brane approaches a\ from below (the (+)-brane can- 
not be located in the southern region a < a\ as discussed 
earlier in subsection B). 

Before closing this subsection, we give the scale factor 
and the scalar field in terms of the cosmic time. In the 
northern region (i.e. a > ax), we find 



($ 4 G 2 - (Br 1 + B2 1 ) 2 ) ($ 4 G 2 - (Sf 1 - B^ 1 ) 2 ) 
/ 2 d$ / d a G 9 a /xx 

+ -w)r (3 ' 46) 



Oi 




(ri 



a/2 



(3.51) 



\a—ai 



This equation can be seen as a first order differential 
equation for the function $(a). Given some initial con- 
dition = $(0*) at some initial point a*, one can in- 
tegrate this differential equation and thus obtain $(a). 
This solution ^(a) can then be substituted into (|3.43[) . 
from which one can extract the evolution of the scale fac- 
tor as a function of the cosmic time r. So far, we have 
not used the second equation (|3.40[) . Since the evolu- 
tion of $ and of a can be determined solely from the two 
other equations, this additional equation plays the role 
of a constraint. This means that the potential cannot be 
arbitrarily chosen but is instead dictated by the geome- 
try. Note that this is a consequence of requiring that the 
bulk is described by a given static solution. Thus, once 
the cosmological evolution is known, the potential for the 
radial scalar field is determined by 



where n the time when the brane reaches i«i . Note that 
B^ 2 > B^ 2 in the northern region (wi < a/p + < 1) 
because R2 < R\ {R+ < Ro for the (+) brane or R < 
i?_ for the (— ) brane). 

The scalar field evolution is then given by 



$ ~ $1 1 + VG'KXn " r) 



V 



1 



2 H \da-) 



$> 2 K(a). 



(3.47) 



which can be expressed in terms of <f> by inverting the 
function $ = $(a). 



(3.52) 



This shows that the scalar field velocity approaches a 
constant and that the ratio $/$, in this limit, is inde- 
pendent of the initial conditions. This analysis is similar 
in the limit a — > a^ on the southern side: one just re- 
places a — a\ by ai — a and B^ 2 — B^ 2 by B^ 2 — B^ 2 ■ 



E. Behavior around the poles 

1. (+)-brane 

In order to analyse the behaviour of the differential 
equation for <J> , Eq. (|3.46p . it is convenient to introduce 
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the function C+(a) denned by 

$ 2 = C| /2 / 1/2 G- 



(3.53) 



If one assumes that C+ is smooth and non-vanishing in 
the limit a — > p+ , one finds rom Eq. (|3.46ll the following 
behaviour: 



1 



c 4 



A 



(3.54) 



where A + is an integration constant and where we have 
introduced the notation 



1 

<?+ 
1 

r+ 



and 



r+ := 



(l-q)V + / 1 
4 \R+ 

(l-q) 2 p + / 1 
4 Vi?+ 



1 



-) 2 



1 



a 3 (l-a 5 ) 



1 -f(P+) ( 1 _ 

32 <?(/?+) V+ g+ ; 20e 2 R+Ro 1 - a 8 

3c (l-a)a 3 (l-a 5 ) 



e 2 £ 2 (5-8a 3 + 3a 8 ) 



(3.56) 



using (|3.1ip in the last equality. 

Note that the condition C + = l/q+ is equivalent to 
the codimension 2 limit 



2tt* V7(a)£+ - Sn 



(3.57) 



where 8n is the deficit angle at the northern pole p = P± 
which can be related to the coefficients cq and c+ [HI, 



cp fry 

C + 27T 



(3.58) 



In order to express the scale factor in terms of the cosmic 
time r, the simplest way is to use Eq. (|3.43[) . which yields 



4C* 4 



4 V+ q + ' 

Then, expanding around C ~ 1 /q + and substituting the 
behaviour Eq. (|3 . 54[) into the above expression, we obtain 



a ~ B+(t+ - r) 



(3.60) 



where 



/3+ = 



S 4 



1-1 + 

^ + (i-^)(-/'(p + ))A + ' 



i/(i-r+) 



4/3 2 



(3.61) 



Thus for r + < 1, we obtain (3 + > 2, namely the velocity 
and acceleration are vanishing at the position of the con- 
ical singularity. For the (+)-brane, we have the condition 
c + > c (R+ < Rq) and thus 



^ 60a 3 (l-«) 2 (l-a 3 )(l-« 5 ) ^ 1 



e z fr + < 



(5 - 8a 3 + 3o^ 2 



< ^ • (3-62) 



One can thus get r + < 1 for sufficiently large e 2 £ 2 . Not- 
ing that 



H 



i-r. 



(i-r+)/2 



p+ - a 



(i+r+)/2 



(3.63) 



we find that, if T + < 1, the brane reaches the northern 
pole in a finite time. Indeed, 



(At) 



p+ 



da 



a, H(0,*(0)) 



(3.64) 



is finite if (1 + T + )/2 < 1 (a* is the initial position of the 
brane). Thus, it is impossible to obtain an ever expand- 
ing brane Universe. 



2. (-)- 



The analysis of the opposite limit a — > p_ = ap + is 
quite similar to the previous case. In the region p_ < a < 
wip + , we always have G(a) < 0, and it is now convenient 
to introduce C_ (a) defined by 



$ 2 = C_(a) 1 / 2 / 1 / 2 (-G)- 1 



(3.65) 



Assuming C_ is non-vanishing and regular in the limit 
a = p-, one can easily find the behaviour 



q- 



A. 



(3.66) 



where A- is an integration constant and 

1 f'(p+a) , 1 1 \ 3 1 1 - a 5 



T_ := 



32 -G(p+a) v r_ q_ ' 
3c (l-a)(l-a 3 ) 
^a<5o7^8^T3)' 



20 e 2 R-R a 5 (l-a 8 ) 
(3.67) 



with g_ and r_ defined as in (|3.55p by replacing R + by 

Note that the condition C_ = l/q~ is equivalent to 
the conical limit 



(3.68) 



where 5g is the deficit angle at the southern pole and we 
have used the relation EH 



c_ 



1 



2tt 



(3.69) 
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The expression of the scale factor in terms of the cosmic 
time is given by 

a-p- ~ B_(r_ - T f- (3.70) 

with 



and we find, as near the northern conical singularity, that 
the brane approache p_ in a finite time if T_ < 1. 

F. Numerical example 

In order to illustrate the cosmological behaviour of the 
branes and to confirm our analytical estimates in the 
various limits discussed above, we have integrated nu- 
merically (|3.46p for the (+)-brane by starting from some 
initial condition = <j>(a»). Our choice of parameters 
is the following: 

a = 0.1, ep+ = l, R+/p+=A.O, R /p+ = 5.0 

(3.72) 

This gives w\ = 0.11696. The initial position and initial 
amplitude are set to be 

M 2 

a., = 0.2p + , $ = 0.01 ^= . (3.73) 

e^/p+ 

From this initial position, we have integrated for decreas- 
ing values of a until reaching w\ as well as for increasing 
values of a up to the northern conical singularity. 

We have plotted the profile of $(a) in Fig. 2. We 
observe that the scalar field vanishes at the conical sin- 
gularity and approaches a constant value, with vanishing 
speed, at w = Wi, which is consistent with our analytical 
solutions. 

In Fig. 3 we have plotted the effective potential and 
total energy densities, respectively \fJV and v7^- Near 
the conical singularity, as shown in Fig. 4, we see that the 
special relation £ = 2V is satisfied, which corresponds to 
the conical singularity condition. 

Our numerical example corresponds to the case where 
the brane reaches the conical singularity in a finite time 
as discussed earlier. 




0.2 0.4 0.6 0.8 1 



FIG. 2: Numerical plot of 3>(a) in the Einstein-Maxwell 
model, with the parameters a = 0.1, ep+ = 1, R+/p+ = 
4.0, Ro/p+ — 5.0. 5>(a) is shown as a function of a (in 

M 2 

units of p+) for an initial amplitude <E>* = 0.01 ^—^= at 
the initial position a, = 0.2p+. 




FIG. 3: Numerical plots of v 7 / V(<&) (blue, dashed curve) 
and \/f S (red, solid curve) are shown as functions of a 
(in the unit of p+) for the solution of Fig. 2. Note that 
the total energy density is vanishing at a = p+wi while 
the potential V r ($) is negative at the same point. 




FIG. 4: Enlargement of Fig. 3 near the northern pole 
(conical singularity) a — p+. As we see, the relation 
E = 2V is satisfied at the conical singularity limit. 



IV. EINSTEIN-MAXWELL-DILATON MODEL 

A. Bulk spacetime in 6D Einstein-Maxwell-dilaton 
model 

In this section, we consider a slightly different model 
based on solutions of the 6D Einstein-Maxwell-dilaton 



theory, which now include a dilaton denoted tp. Our bulk 
action is based on the bosonic part of the Nishino-Sezgin 
supergravity [ill ]: 

S B = J d 6 x^(±R - l -{d A ^f - l e -*F AB F AB 
, (4.1) 
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where the other fields can be consistently set to zero. It 
is possible to find exact solutions, given by [ll[ 



ds 2 = {2p + w) r q flu dx^dx 1 ' 
p + / Adw 2 
+ 2g^\(l-a) 2 f(w) 



+ f(w)c 2 d8 2 ), (4.2) 



F,„fl — — 



V2c c 



where 



f(v> 



(1 - a)g w 3 

_ 2(1 - w 2 )(w 2 - a 2 ) 
w 3 (l — a) 2 



<p = -hx(2p+w), (4.3) 



(4.4) 



The geometry possesses two conical singularities at w = 
a and w = 1, respectively. As in the non-dilatonic case, 
we regularize the corresponding codimension 2 branes by 
introducing two codimension- 1 branes at the respective 
positions w- > a and w + < 1, delimiting two regular 
caps. The northern cap (w+ < w < 1) and the southern 
cap (a < w < w_), can be described by solutions similar 
to that of the main bulk, i.e. of the form (14.2|) - (|4.3p . but 
with parameters g± and c± that differ from go and Co- 
The regularity of the caps imposes [12] 



1 



1 + a 1 + a 

while the continuity of the (6, 9) component requires 



(4.5) 



2 9 2 + 2g 2 _ 2g 2 



(4.6) 



The solutions for the bulk gauge field around the north 
and south poles are given by 



.4 



1 — a ur 1— a ur or 



The difference should be a pure gauge term, which im- 
plies the Dirac quantization condition 



N = 



e£(l + a) 



a 



(4.8) 



B. Junction conditions and brane dynamics 



We describe the brane matter by an action of the form 



Sh.+ = — 



$ 2 ± (d a a± - eA a 



(4.9) 



It is similar to the non-dilatonic case with two differ- 
ences: the potential V($±,(p) now depends on the dila- 
ton, and we have introduced a coupling £ (ip) , assumed to 
be strictly positive, between the kinetic term of the scalar 
field and the bulk dilaton. The energy density E and 



the pressures along the external and internal directions, 
respectively p and p s , are defined like in (|3.15[) - (|3.17[) . 
with the only difference that the potential now depends 
on the dilaton if, in addition to the scalar field The 
comparison of the metric (|4.2[) with Eq. (|2.ip yields the 
identification 

2p + 



A 2 



2p + w , B l 



P+ 2 f( \ 
— 2 C J/W 



2gj 



g 2 (l-a) 2 f(w) 



(4.10) 



and the scale factor of the brane is, as before, a(r) = 
A(wt(r)). As in the previous section, we also define the 
rescaled metric component 



(2 P+ ) 1 / 2 V7M 



a 



(4.11) 



The junction conditions for the metric are exactly the 
same as in the previous section, and are thus given by 
p.l5H3.l7)l . They imply in particular the generalized 
Fricdmann equation 

a 2 p 2 



167r 2 Pp+f(3 + i 7 



7^2 fp + (3 + 



2/' 



Bo 2 ) 2 



B n 



P 2 



2a 2 



(4.12) 

'y/fY.. The Friedmann equation is 



where p :— 2n£^/p~\ 
very similar to the one obtained in the Einstein-Maxwell 
case and one would find, by repeating the same procedure 
as in the previous section, that the low-energy behaviour 
differs from standard cosmology. 

The Maxwell junction condition, however, is modified 
because the dilaton is explicitly coupled to the electro- 
magnetic field. It now reads 



"^(a) 



X 



-e£,(<p)<5> 2 (dga-eA e 



where we have defined 



F w g — — 



16a£ pi 



1 



(4.13) 



(4.14) 



Moreover, we have an additional junction condition, that 
of the dilaton, which is given by 



X 



<p(a) 



2 
a 



eA f 



(4.15) 



As in the Einstein-Maxwell model, the bulk is divided in 
two by a limit that cannot be crossed by the branes. In 
the present case, the position of this limit is given by 



Wl 



(4.16) 
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Similarly, the total energy density in the brane is negative 
when 



a" jl + a 2 

wi < - — < w 2 
2p + 



(4.17) 



and positive otherwise. As in the case of the Einstein- 
Maxwell model, the (+)-brane cannnot be in the region 
(2p_) 1/2 < a < (2,0+wi) 1 / 2 , where p- := p+a. The 
brane can also cross W2 simply with a change of the sign 
of the energy density. 

Following the same method as in the Einstein-Maxwell 
case, one can rewrite the junction conditions as 



B~ 2 + a 2 - ^B- 2 + a 2 = £(^)$ 2 G(a) , (4.18) 
^ (^M* 2 + V($, <pj) = ^)K{a) , (4.19) 



B- 2 + a? xB- 2 +a? 



a^ 2 



(4.20) 



where 
K(a) 
G(a) 



e 2 (l-a) 2 t ,_ s _ 6 fdJ , 7\(d a f 1 



64o! 2 (2p + ) 3 '' K ~'~ V 2/ ' aJ\2f a 
/(«)«•(*■> 1 



(4.21 



32a 2 (2p + ) 3 '' w V 2f a 
This implies in particular 



One can also note that the topological conditions are 
the same as in the static case [8] (see Appendix A-2 for 
details). 

To obtain the cosmological evolution, one can proceed 
exactly as in the previous section. In particular, the 
square of the scale factor velocity is given by 



" 2 = 4^2^)2 (<f 4 G 2 £M 2 - (V + B 2 -1 ) 2 ) ~ (Sf 1 - B a -1 ) 2 ) := *(«)) 2 . (4-23) 



and the differential equation for $ is now 



f^ 2 = 4^G 3 ^) 3 / 2 / ^G tffo) 

Ida/ ($ 4 G 2 ^) 2 -(^ 1 +5 2 -V)(* 4 G 2 e(^) 2 -(S 1 " 1 -S 2 - 1 ) 2 ) V$da VG a K y/J ' J J y ' ' 



The junction condition (|4.19[) implies 



scalar field Eq. (|4.24|) becomes in the limit <I> 4 G 2 £ 2 — > 



"-eMC-Ms)'**"*")- <4 ' 25 > 

which provides an expression for the potential y as a 
function of a. The potential depends on both the scalar 
field $ and the dilaton <p, and it is possible to disentangle 
the dependence on $ and if by using the dilaton junc- 
tion condition Eq. (|4.15|) with Eq. (|4. 18|) . which yields 
the following expression for the partial derivative of the 
potential with respect to the dilaton: 



dV 

dtp 



-*(a) 2 G(u)(l(( v («)) - \(' («>(«))(£ - ;)) 



. da / 



W(o,#(o))" 



(4.26) 



C. Behaviour near the "barrier" 

We now analyze the behavior near the critical point 
a — > ai := (2/J4.1U1) 1 / 2 . The differential equation of the 



d$\2 

5 / 

8$ 6 G 3 £ 3 



b: 



B- 2 ) 



rid* lfd a G (<*\\\ 

2 U^ + 2{— + daln {7?))r 

(4.27) 



We expand G(a) — G'(ai)(a—ai) andi^(a) = K'{a\)(a— 
oi), with G'(ai) > and K'(a\) < 0. Then, as- 
suming that £ tends smoothly towards a finite value 
£i := £{tp(a>i)), we find the following approximate be- 
haviour 



$ ~ $o + $i(a - ai) 2 
The potential is given by 

V~-^ZZ(V)G' 



$2= e( V ) a $gG' 3 



^2 



(B 2 "-r) 



< 0, 



(4.28) 

a— ai 



(4.29) 



where we have used Eq. (|4.28[) . Thus, the potential is 
finite and negative in the limit a—*ai. 
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We now derive the expression for the scale factor in 
terms of the brane proper time. In the case a — + , we 
obtain 



- (2p + w 1 ) 1 / 2 = 



B 2 2 - B l 



a=(2p + to 1 ) 1 /2 



(n-r) 



where 



i-r. 



'?+(£- £)A+ -/'((2p+) 1/2 ) 



(4.30) 



4/3i 



(2P+) 



i/(i-r+) 
(4.38) 



where ri represents the time when the brane reaches a%. 
Note that B^ 2 > B^ 2 for cii < a < (2p + y/ 2 because 
92 > .91 (<7+ > go)- The scalar field configuration is given 
by 



$ = $o(l + v / G 7 ?M 



a=(2 P+ w 1 ) 1 / 2 



n - r) . (4.31) 



Thus the scalar field has a constant velocity. The brane 
terminates its motion at the point. The analysis a — > 
can be done in the similar manner: one just replaces 
a — a,\ by cii — a and -B^ -2 — B^ 2 by 2 — B^ 2 ■ 



D. Behaviour around the poles 

In the limit a — > (2/5+) 1 / 2 , it is useful to rewrite the 
scalar field as 



$ 2 = C + (a)^^ G -i. 



(4.32) 



If C+(a) is a smooth and non- vanishing function near 
a = (2P+) 1 / 2 , an approximate solution of the differential 
equation is given by 



1 

9+ 



C + ^A + {{2p + Y' 2 ~a)^ 



(4.33) 



where 



1 

9+ 
1 

r+ 



(1 - «) 2 P+ 
2 

(1 - a) 2 p+ 



(.9+ - 9oY 

{g+ + go)' 



(4.34) 



and 



r+ = 



1 -/'((2p + )V») , i 



32 a 2 £G 



(2P+) 1 



1 



a cq 



r + 9 + / e 1 + a 



(4.35) 



Note that the condition C+ = I/9+ is equivalent to the 
conical limit 



27r v /pPv / 7R£ + 



= <5 W , 



(4.36) 



where Sn represents the deficit angle at the northern 
pole. The expression for the scale factor in terms of r 
is given by 



(2p+) 



1/2 



B + (r + -rf\ 



(4.37) 



For the (+)-brane, we have the condition c + > cq 
(.9+ > go) an d a discussion similar to that of the Einstein- 
Maxwell model tells us that it is possible, with a suffi- 
ciently large charge for the brane, to have T + < 1, in 
which case the brane reaches the northern pole in a fi- 
nite time. Thus, one does find an ever expanding brane 
Universe. 

The results for the opposite limit a — > (2p_) 1 / 2 are 
quite similar to the above results. Using the decomposi- 
tion, 



$2 



C_(o) 



1/2 



1/2 



one can identify the behaviour 



(4.39) 



9- 



A- 



(a-(2p_) 1 /2) r - (4.40) 



where 



1 f((2p-) 1/2 ) 
32 -a 2 £G| (2p )1/2 



1 



1 

9- 



1 



Co 



e 2 ^ 2 x + a 

(4.41) 



and 9_ and r_ are defined as in (|4.34[) by replacing 
with g_ . 

Note that the condition C_ = I/9- corresponds to 



2ttVpPV7W£- 



(4.42) 



where 8s denotes the deficit angle at the southern pole. 
In terms of the cosmic proper time t, we find 



a - (2p_) 1/2 = B_(t_ - t 



(4.43) 



where 
0- 
B_ 



i-r_ 



4/3 2 



(2p- 



i/(i-r_) 

(4.44) 



For the (— )-brane, c_ > Co (.9- > go) aud one can have 
r_ < 1, in which case the brane reaches the south pole 
in a finite time. 
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V. DISCUSSIONS 

In this work, we have investigated the dynamics of the 
regularized branes in a 6D bulk. We considered models 
based on either the Einstein-Maxwell theory or Einstein- 
Maxwell-dilaton theory (more precisely a bosonic part of 
Nishino-Sezgin supergravity) . In both models, the origi- 
nal systems are composed of a warped bulk bounded by 
codimcnsion 2 branes located at the poles. In order to 
introduce non-trivial matter on the branes we have reg- 
ularized the codimension 2 branes, located at the poles, 
by replacing them with ring-like codimension 1 branes. 

The brane matter is composed of a complex scalar field 
coupled to the bulk U(l) gauge field (and to the bulk 
dilaton in the dilatonic model). Such matter was used 
for a static (fixed) brane, where the radial mode of the 
complex scalar field was assumed to be stabilized at a 
local minimum of its potential and thus frozen at a fixed 
vacuum expectation value 0. For moving branes, one 
must allow for a time dependence of the radial mode and 
this radial mode controls the dynamics of the brane in 
the bulk. 

We have studied the dynamics of such branes in the two 
types of models mentioned above. In both models, branes 
exhibit very similar behaviors and their essential features 
can be summarized as follows. First, there is a critical 
radius, translated into a critical scale factor a\ given by 
W\p+ in the Einstein-Maxwell model (see Eq. (|3.37p ) 
and y/'2p + wi in the Einstein-Maxwell-dilaton model (see 
Eq. (|4.16p ). which cannot be crossed by a brane: an 
expanding brane with a < a\ at some time cannot expand 
beyond ax', conversely, a contracting brane with a > a\ 
cannot shrink below a±. Investigating the motion of a 
brane close to this critical point, we have found that a 
brane reaches the point a±, within a finite time, with a 
divergent velocity (d cx |t — ti| -1 / 2 , where r is the proper 
time on the brane). Therefore, the Hubble parameter of 
the brane diverges at a = a\. Second, near the poles, 
we have found behaviours where the brane reaches the 
conical singularity within a finite time with vanishing 
velocity and acceleration. This occurs if the brane charge 
is sufficienly large. 

For the Einstein-Maxwell model, we have solved the 
brane motion numerically and confirmed our analytic re- 
sults in both limits. In both types of models, the size 
of the internal dimension goes to zero as the brane ap- 
proaches a conical singularity or to a constant as the 
brane approaches the critical value a\. These conclu- 
sions were obtained by assuming a complex scalar field 
as brane matter and it would be worthwhile investigating 
the cosmological behaviour for other kinds of matter. 

More generally, i.e. independently of the specific brane 
matter content, we have seen that the Friedmann equa- 
tions Eq. (|3.24j) and (|4.12j) . in the Einstein-Maxwell and 
Einstein-Maxwell-dilaton models, respectively, do not be- 
have like standard cosmology in the low energy limit, be- 
cause the effective gravitational coupling, Eq. (|3.28|) is 
strongly time-dependent (and even negative in the vicin- 



ity of the pole). Thus, our results show that, in con- 
trast with the 5D Randall-Sundrum cosmology in the 
low-energy limit, the dynamics of the brane in the present 
setup is not compatible with the usual cosmological evo- 
lution and modifications are necessary to obtain a sen- 
sible cosmology. One possible modification would be to 
consider time-dependent bulk solutions [2l[ instead of the 
static geometries considered here. 2 . We will report these 
investigations in future publications. 



VI. ADDENDUM 

While this paper was being completed, we became 
aware of a related work (23|, which studies the cosmo- 
logical behaviour in the same type of models as those 
discussed here. As far as we can see, our results agree 
with those of [23| . The two works are also complementary 
in many respects. 
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APPENDIX A: BRANE EQUATIONS AND 
QUANTUM NUMBERS 

1. Einstein-Maxwell model 

The dynamics of the radial scalar field on the brane is 
given by 

* + Ki + c^) * + y * + h ( dea ~ eAe ) 2$ 

= 0. (Al) 

The equation of motion of <& field is consistent with the 
energy momentum conservation law on the brane. 



fi + ^( s+P 

ith 



C'A 
~C~ 



AR 



Tabu n 



S + PB ) = 



= A [Fa 



t A I 

Tabu n 



X 



= AC- 2 e<$> 2 (d e o-eAe)F A 6 



,(A2) 



(A3) 



2 Low energy behaviour of brane gravity in the Einstein-Maxwell 
model was discussed recently in Ref. 1221 
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where, in the final step, we have used the Maxwell junc- 
tion condition. 

The equation of motion for the Goldstone modes on 
the branes is given by 



d a (y=-q<S> 2 q ab (d b a-eA b )) = 0, 



(A4) 



which implies d 2 a± = 0. The equation of motion for 
the axial mode a can be solved as a± = n±9. From Eq. 
p,12p . (|3.29|) . (I3.31|) . one finds that the brane quantum 
numbers n± are given by 



n+ 



N 5 - 8a 3 + 3a 8 
Y4(l-a 3 )(l-a 5 ) ' 

N 5-8a- 3 + 3a~ 8 
~Y4(l-a- 3 )(l-cr 5 ) 



(A5) 



where we have also used Eq. ()3.13p . Thus, between quan- 
tum numbers, we have the following relations: 



N. 



(A6) 



Note that the same relations as Eq. (|A5[) are obtained 
in the case of the static brane. Due to the redundancy 
of the junction condition, they appear as a constraint 
between the model parameters. 



2. Einstein-Maxwell-dilaton model 

The dynamics of the radial scalar field on the brane is 
jiven by 



+ + $ + ?<j& + V* + -^ (d e a - eA ( 



The equation of motion of <& field is consistent with the 
energy momentum conservation law on the brane: 



where 



C'A 



S + PB ) = 



n- A I 

Tabu n 



.{AS) 



Tabu n 



a„b 



= A 



X 



= -e64C" 2 $ 2 (d e a - eA e )e- <p F AB + <p, A V, v A 

+ ^{<p)<pM( - & + ® 2 C- 2 (dev - eAe)) , 

(A9) 



and at the final step we used the Maxwell and dilaton 
junction conditions. 

The equation of motion for the Goldstone modes on 
the branes is again given by 



d a {^qav)® 2 q ab (d b v-eA b )} =0, (A10) 
which implies dg(T — 0. The solution is also given by 
a = n±6. The quantum numbers n± are found to be 
given by 



n± 



= 0. 



4> 
(A7) 



N 



N . 



(All) 
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